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We examine how nuclear masses are related to the density dependence of the symmetry energy. 
Using a macroscopic nuclear model we calculate nuclear masses in a way dependent on the equation 
of state of asymmetric nuclear matter. We find by comparison with empirical two-proton separation 
energies that a smaller symmetry energy at subnuclear densities, corresponding to a larger density 
symmetry coefficient L, is favored. This tendency, which is clearly seen for nuclei that are neutron- 
rich, nondeformed, and light, can be understood from the property of the surface symmetry energy 
in a compressible liquid-drop picture. 

PACS numbers: 21.65.Ef, 21.10.Dr 



INTRODUCTION 



o 

(N 
u , 

< 

Saturation of the density and binding energy is a fundamental property of atomic nuclei. Due to this property, 
the nuclear masses can be described well using a liquid-drop approach. Conventionally, in this approach, the nuclear 
binding energy Eb is written as function of mass number A and charge number Z (or neutron number N) with the 
Weizsacker-Bethe mass formula [l| 

O ■ — Eb = Eyol + Esym + -Esurf + -Ecoul, (1) 

' where E'voi = flvoi^ is the volume energy, £^sym = asym[(-^ — Z)/A]'^A is the symmetry energy, E'surf = Osurf^^^^ is the 
surface energy, and i?coui = acouiZ^/A^^^ is the Coulomb energy. The sum Eyo\ + Egym corresponds to the saturation 

[ energy of uniform nuclear matter. Since the matter in a nucleus constitutes a strongly interacting system, it remains 

^ . a challenging theoretical problem to understand the nuclear matter equation of state (EOS) through microscopic 

■ calculations that utilize a model of the nuclear force duly incorporating low-energy two-nucleon scattering data and 

CN| , properties of very low-mass nuclei 0- Furthermore, it is not straightforward to empirically clarify the EOS, although 

CO ■ constraints on the EOS obtained from nuclear masses and radii (e.g., Refs. [2-B)i_pbservables in heavy- ion collision 

I experiments performed at intermediate and relativistic energies (e.g., Refs. [3, @), the isoscalar giant monopole 

ps) ■ resonance in nuclei (e.g., Ref. 0), and even X-ray observations of isolated neutron stars and quiescent low-mass 

' X-ray binaries [Tlj do exist. In this work we will consider such constraints provided by masses of unstable nuclei. 

The energy density of bulk nuclear matter is a function of nucleon density n and proton fraction x, which are related 

I ' to the neutron and proton number densities n„ and rip as n„ = n{l — x) and nj, = nx. We can generally express the 

^ . energy per nucleon near the saturation point of symmetric nuclear matter as |12l | 
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Here wq, ?^o and Ko are the saturation energy, the saturation density and the incompressibility of symmetric nuclear 
matter, and a = 1 — 2x is the neutron excess. L and Sq are associated with the density dependent symmetry energy 
coefficient S{n): Sq is the symmetry energy coefficient at n — no, and L = 3no{dS / dn)n=na is the symmetry energy 
density derivative coefficient (hereafter referred to as the "density symmetry coefficient"). As the neutron excess 
increases from zero, the saturation point moves in the density versus energy plane. This movement is determined 
mainly by the parameters L and Sq. Up to second order in a, the saturation energy Wg and density ris are given by 

Ws ^ Wo + Soa^ (3) 

and 

Us^uo -r^a . (4) 

The slope, y, of the saturation line near a — Q {x ^ 1/2) is thus expressed as 
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In our earlier investigations f6| we explored a systematic way of extracting L and S'o from empirical masses and 
radii of nuclei, together with the parameters, no, Wq and Kq, characterizing the saturation of symmetric nuclear 
matter. We first gave an expression for the energy of uniform nuclear matter, which reduces to the phenomenological 
form ^ in the simultaneous limit n — >■ no and a — t- (x — >■ 1/2). Using this energy expression within a simplified 
version of the extended Thomas-Fermi approximation, which permits us to determine the macroscopic features of 
the nuclear ground state, we calculated charges, charge radii and masses of /3-stable nuclei for fixed A. Comparing 
these calculations with empirical values allows us to derive the optimal parameter set for fixed values of the slope y 
and the incompressibility Kq. We thus found a strong correlation between L and S'o. The next step was to calculate 
the root-mean-square (rms) charge and matter radii of more neutron-rich nuclei that are expected to be produced in 
radioactive ion beam facilities. The results suggest that the density symmetry coefficient L may be constrained by 
possible systematic data for the matter radii in a manner that is nearly independent of Kq. 

The reason that we concentrated on radii of unstable nuclei rather than their masses in Ref. [61 was that we 
originally considered the surface and electrostatic properties as well as shell and pairing effects, which play a role in 
nuclear masses, to obscure the derivation of the EOS parameters. Nevertheless, in the present analysis, we perform a 
systematic calculation of nuclear masses from the same framework used for our previous calculations of nuclear radii. 
Since experimental mass data have been accumulated even for unstable nuclei, we can compare the calculations with 
the existing data, which is a great advantage over the case of nuclear radii. Comparison of the two-proton separation 
energy {S2p{Z,N) — Eb{Z,N) — Eb{Z — 2, A^)) implies that a larger i, corresponding to smaller symmetry energy 
at subnuclear densities, is favored. This tendency is most clearly elucidated from the masses of light, nondeformed, 
and neutron-rich nuclei. This can be understood from a compressible liquid-drop picture of nuclei in terms of the L 
dependence of the surface symmetry energy, which has some relevance to the neutron skin thickness [Tsl . [l^ . 

In Sec. II we summarize a macroscopic model of nuclei used here. Calculations of nuclear masses and the comparison 
with empirical data are illustrated in Sec. III. Our conclusions are presented in Sec. IV. 



II. MACROSCOPIC NUCLEAR MODEL 

In this section, we summarize a macroscopic model of nuclei @ , which was constructed in such a way as to reproduce 
the known global properties of stable nuclei and can be used for describing the masses and radii of unstable nuclei in 
a manner that is dependent on the EOS of nuclear matter. 

The bulk energy per nucleon is an essential ingredient of the macroscopic nuclear model. We set this energy as 

lUm„n 

+ {1 — a^)vs{n)/n + a^Vn{n)/n, (6) 
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are the potential energy densities for symmetric nuclear matter and pure neutron matter, n„ and n^ are the neutron 
and proton number densities, n = n„ -I- np, a = (rin — np)/n is the neutron excess, and m„ is the neutron mass. 
Expressions ([6])-(|8]) can well reproduce the microscopic calculations of symmetric nuclear matter and pure neutron 
matter by Friedman and Pandharipande in the variational method. In this method, the isospin dependence of 
asymmetric matter EOS is shown to be well approximated by Eq. © [l6j . (Replacement of the proton mass nip by 
TO„ in the proton kinetic energy would make only a negligible difference.) For the later purpose of roughly describing 
the nucleon distribution in a nucleus, we incorporate into the potential energy densities ([7]) and ([8]) a low density 
behavior oc n^ as expected from a contact two-nucleon interaction. A set of expressions ([S])-® is one of the simplest 
that reduces to the usual form Q in the limit of n — > no and a — > 0. In fact, the parameters ai, • • • , 63 are related to 
no, wq, Kq, Sq, and L as 
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We determine the parameters oi, • ■ • , 63 in such a way that the charge number, charge radius, and mass of stable 
nuclei calculated in a macroscopic nuclear model constructed in Ref. [6| are consistent with the empirical data. In the 
course of this determination, we fix 63, which controls the EOS of matter for large neutron excess and high density, 
at 1.58632 fm'^. This value was obtained by one of the authors in such a way as to reproduce the neutron matter 
energy of Friedman and Pandharipande [15j. Change in this parameter would make no significant difference in the 
determination of the other parameters and the final results for nuclear masses. 

We describe macroscopic nuclear properties in a way dependent on the EOS parameters ai,---,63 by using a 
Thomas- Fermi model Q. The essential point of this model is to write down the total energy of a nucleus of mass 
number A and charge number Z as a function of the density distributions 7i„(r) and 7ip(r) in the form 



E = Eb + E„+Ec + Nm-ac^ 



ZrripC 



where 



is the bulk energy, 



is the gradient energy with adjustable constant Fq, 
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is the Coulomb energy, and N — A — Z is the neutron number. This functional allows us to connect the EOS and the 
nuclear binding energy through the bulk energy part Et- For simplicity we use the following parametrization for the 
nucleon distributions ni{r) [i = n,p): 



ni{r) = 



Ri 



0, 



r <R, 



r > R, 



(18) 



This parametrization allows for the central density, half-density radius, and surface diffuseness for neutrons and 
protons separately. 

In order to construct the nuclear model in such a way as to reproduce empirical masses and radii of stable nuclei, 
we first extremized the binding energy with respect to the particle distributions for fixed mass number, five EOS 
parameters oi, • • ■ , &2j and gradient coefficient Fq. Next, for various sets of the incompressibility Kq and the density 
symmetry coefficient L, we obtained the remaining three EOS parameters 710, wq, and So and the gradient coefficient 
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FIG. 1: (Color online) The sets of {L,Ko) (crosses) consistent with the mass and radius data for stable nuclei. The thin 
lines are lines of constant y. The labels A-I denote the sets for which we performed detailed calculations of the ground state 
properties of inhomogeneous nuclear matter at subnuclear densities . In this paper we often focus on the EOS models with 
the parameter sets C and G. For comparison, the values calculated from two mean-field models [TMl (square) and SIII (dot)], 
which are known to be extreme cases are plotted. The plot shows that our sets of {L,Ko) effectively cover such extreme 
cases and constraints on L and Ko from other observables [l9l ]. 




Fq by fitting the calculated optimal values of charge number, mass excess, rms charge radius to empirical data for stable 
nuclei on the smoothed f3 stability line [l^ and by using Eqs. (Pjl- P^ . In the range of the parameters < L < 160 
MeV and 180 MeV < Kq < 360 MeV, as long as y < —200 MeV fm'^, we obtained a reasonable fitting to such data 
(see Fig. 1). As a result of this fitting, the parameters no, wq, So, and Fo are constrained as no = 0.14-0.17 fm~^. 
Wo = — 16 ± 1 MeV, 5*0 = 25-40 MeV, and Fq = 66 ± 6 MeV fm^. We remark that a negative L is inconsistent with 
the fact that the size oi A = 17, 20, 31 isobars deduced from the experimental values of the interaction cross section 
tends to increase with neutron/proton excess [l^. This inconsistency can be seen from Eq. ^ which shows that the 
saturation density Ug increases (and hence the isobar size decreases) with neutron/proton excess for a negative L. We 
also note that the fitting gives rise to a relation nearly independent of Ko, 

Sa^B + CL, (19) 

where i? w 28 MeV and C « 0.075. As we shall see, this linear relation plays a part in the L dependence of calculated 
masses at large neutron excess. 

We remark that in the range of the EOS parameters L and Kq shown in Fig. 1, the calculations agree well with 



a more extended data set of nuclear masses for A > 2 [2l| and charge radii for A > 50 [22]. The rms deviations of 



the calculated masses from the measured values are ~ 3-5 MeV, which are comparable with the deviations obtained 
from a Weizsacker-Bethe formula, while the rms deviations of the calculated charge radii from the measured values 
are about 0.06 fm, which are comparable with the deviations obtained from the A^^'^ law. As we shall see, detailed 
comparison with empirical masses in terms of two-proton separation energy allows us to tell which of the EOS models 
with large L and with small L are more favored. 



III. NUCLEAR MASSES 



We now proceed to evaluate nuclear masses from various EOS models with the parameter set {L,Ko) consistent 
with the mass and radius data for stable nuclei by minimizing Eq. (ITi)) for fixed N and Z. The results are then 
compared with the existing data in terms of the two-proton separation energy. We finally discuss the L dependence 
of the calculated masses within the framework of a compressible liquid-drop picture of nuclei. 

We begin by illustrating the calculated and experimental values of the two- proton separation energy S2p for O, Mg, 
Ca, Ni, Sn, and Pb isotopes. The two-proton separation energy is useful partly because the even-odd staggering is 
essentially cancelled out and partly because the isotope dependence except for shell gaps is smooth according to the 
Yamada-Matumoto systematics [23]. In fact, as can be seen from Fig. 2 in which the two-proton separation energy 
minus the one calculated from a Weizsacker-Bethe mass formula ([T]) with the coefficients Ovoi = —15.5391 MeV, 
flsym = 22.7739 MeV, a^urf = 16.9666 MeV, and acoui = 0.703893 MeV [H, i.e., S'zp - S^^, is plotted, the empirical 
values behave monotonously with neutron excess except for the vicinity oi N — Z and neutron magic numbers and 
for the deformed region. Note that this monotonous behavior is clearer for O, Mg, and Ca than for Ni, Sn, and Pb. 

We remark that the macroscopic nuclear model summarized in the previous section was originally used for systematic 
calculations of charge and matter radii of nuclei of A ^ 50 Here, we performed systematic calculations of masses 
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FIG. 2: (Color online) The two-proton separation energy, having the one calculated from a Weizsacker-Bethe mass formula [23] 
subtracted out, for O, Mg, Ca, Ni, Sn, and Pb isotopes. The empirical values [2ll ]. the calculated values from the EOS models 
G and C, and the values obtained from a contemporary mass formula [1^ are plotted in each panel. 



of nuclei including O, Mg, and Ca isotopes. Applicability to lighter nuclei can be confirmed from Fig. 2, which shows 
that differences between the calculated and empirical values of S2p for stable nuclei are limited within 3 MeV. 

The calculated values of S2p, neither including deformation, the Wigner term, nor shell corrections, show a smooth 
dependence on neutron excess in a way different between the EOS models G (L = 5.7 MeV) and C {L = 146 MeV). 
(The Weizsacker-Bethe formula, which is based on the incompressible liquid-drop model for nuclei, is independent of 
the parameters L and Kq characterizing the density dependence of the EOS.) In order to examine how the calculated 
values of S2p depend on L and Kq, we plot in Fig. 3 the results obtained for ^^Ni and ^^O from the EOS models 
with various values of L and Kq shown in Fig. 1. The results decrease with L almost linearly, while they are nearly 
independent of Kq. Accordingly, for various values of L one can predict where 5*2^ — is located in Fig. 2 by 
interpolation or extrapolation of the results from the EOS models G and C. As far as the slope of 5*2^ — is 
concerned, the result with the EOS model C seems to be more consistent with the empirical behavior, particularly 
for light nuclei. 

We also notice from Fig. 2 that there is a roughly uniform offset between the empirical values and the values 
calculated from the EOS model C at > Z. This offset comes mainly from the proton shell gaps. In fact, it is 
effectively cancelled out in a region away from Z = N and neutron magic numbers by taking a difference, S2p{Z, N) — 
S2p{Z,N — 2) = ASVnp, as can be seen from Fig. 4. In contrast, the empirical behavior of this difference shows a 
N — Z and neutron shell structure in a more exaggerated manner. As discussed in Ref. [2^, the smoothed behavior 
of SVnp is related to the coefficients agym and assym affixed to A[{N — Z)/A]'^ and A^^^KN — Z)/A]'^ in the mass 
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FIG. 3: (Color online) The two-proton separation energy calculated for Ni and O as a function of L 
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FIG. 4: (Color online) SVnp for O, Mg, Ca, Ni, Sn, and Pb isotopes. In each panel, the empirical values f2l)|, the values from a 
contemporary mass formula [2^ and a Weizsacker-Bethe mass formula [S^], and the calculated values from the EOS models C 
and G are plotted. 



formula ([T]) (with the surface symmetry term added) by 

SVnp ~ 2{asym + assymA-l/3)/^_ (20) 

Our results for SVnp suggest that the L dependence lies in the parameters Usym a-nd assym- In fact, this suggestion 
can be understood within the framework of a compressible liquid-drop model in which nuclei in equilibrium are allowed 
to have a density different from the saturation density uq of symmetric nuclear matter. If one ignores Coulomb and 



7 



-22 
-24 
-26 
-28 
-30 
-32 
-34 
-36 



78 



Ni 



J ±. 



50 



100 
L (MeV) 



150 




100 

L (MeV) 



FIG. 5: (Color online) The mass excess calculated for "*Ni and ^^O as a function of L. 



surface corrections, the equilibrium density and energy per nucleon of a liquid-drop correspond to rig and Ws-, which 
are given by Eqs. ^ and ^ for nearly symmetric nuclear matter. Note that the parameter agym is characterized by 
the symmetry energy coefficient Sq, which in turn is related to L by Eq. (1191) . Since one obtains a larger agym for 
larger L, the effect of agym tends to increase SVnp with L, which is in the wrong direction. This effect, therefore, has 
to be dominated by the effect of agsym- This is consistent with the fact that the L dependence of the calculated S2p is 
clearer for lighter nuclei and also with the tendency that the calculated mass decreases with L for neutron-rich nuclei 
(see Fig. 5). 

The effect of agsym can be understood by considering the density-dependent surface tension [isj . 
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where riin and ain are the density and neutron excess inside a liquid-drop, cfq — o'(no, 0), Cgym is the surface symmetry 
energy coefhcient, and x = inQ/ao)d(T /dnin\nin=no-aia=Q- By substituting Hs into rijn, one obtains 



(T(ns, ain) = (Jo 
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Thus, assym^^^'^ behaves as —AiraoR^ (^Csym + ; with the liquid-drop radius R. The condition x > 0, which is 
suggested by various nuclear models including the macroscopic one utilized here [27| . is desirable for understanding 
of the L dependence of SVnp. We remark that this condition is suggested by various mean- field models because the 
calculated neutron skin thickness, which is basically proportional to — Ossym, increases with L [l3 |. 



IV. CONCLUSIONS 



We have analyzed the influence of the density dependence of the symmetry energy on nuclear masses by using a 
macroscopic nuclear model that depends explicitly on the EOS of nuclear matter. We find that the L dependence of 
the calculated masses comes mainly from the surface property through the density and neutron excess dependence of 
the surface tension. 

In making reasonable estimates of L from empirical masses, we find that empirical two-proton separation energies 
for neutron-rich, light nuclei are useful, implying a larger L value. However, we still have two caveats. First, Eq. 
(im suggests that uncertainties in xl^o affect the prediction of the value L. Second, a smooth isotope dependence 
of empirical mass data is hard to derive. Consequently, such estimates of L are closely connected to determination of 
the parameters x s-^id characterizing the density dependence of the surface and bulk energy and to understanding 
of the discrete behavior of and SVnp [H, HI, Hll • 
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